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ABSTRACT 

During our Herschel Lensing Survey (HLS) of massive galaxy clusters, we have discovered an exceptionally bright source behind the z=0.22 
cluster Abell 773, which appears to be a strongly lensed submillimeter galaxy (SMG) at z=5.2429. This source is unusual compared to most other 
lensed sources discovered by Herschel so far, because of its higher submm flux (~200mJy at 500yum) and its high redshift. The dominant lens is a 
foreground z=0.63 galaxy, not the cluster itself. The source has a far-infrared (FIR) luminosity of L FIR = 1.1 10 14 /yu L G , where p. is the magnification 
factor, likely ~ 11. We report here the redshift identification through CO lines with the IRAM-30m, and the analysis of the gas excitation, based 
on CO(7-6), CO(6-5), CO(5-4) detected at IRAM and the CO(2-l) at the EVLA. All lines decompose into a wide and strong red component, and 
a narrower and weaker blue component, 540km s apart. Assuming the ultraluminous galaxy (ULIRG) CO-to-Hi conversion ratio, the H2 mass 
is 5.8 lO'Vy" Mo, of which one third is in a cool component. From the Cl( 3 />2 _3j Pi) line we derive a Cl/H 2 number abundance of 6 10~ 5 similar 
to that in other ULIRGs. The H2O,,(2,0,2- 1,1,1) line is strong only in the red velocity component, with an intensity ratio I(H20)/I(CO) ~ 0.5, 
suggesting a strong local FIR radiation field, possibly from an active nucleus (AGN) component. We detect the [NII]205/jm line for the first time 
at high-z. It shows comparable blue and red components, with a strikingly broad blue one, suggesting strong ionized gas flows. 
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1. Introduction 

At high redshift, z~l-3, dusty luminous infrared galaxies 
are thought to dominate the history of cosmic star forma- 
tion (Dole etal. 2006; Wardlow etal. 2011). At these epochs, 
star formation was occuring mainly in massive galaxies, a 
tendency dubbed downsizing in galaxy formation scenarios 
(e.g. Cowie etal. 1996; Heavens etal. 2004; Perez-Gonzalez 
et al. 2008; Magnelli et al. 2009, 201 1). These massive starbursts 
might correspond to the formation phase of the luminous el- 
lipticals seen in high-density regions today (Genzel etal. 2003; 
Swinbank etal. 2006). They are therefore essential for under- 
standing the main processes of galaxy formation and evo- 
lution. In particular, atomic and molecular lines are crucial 
clues in these systems for tackling star formation and AGN 
(Active Galaxy Nucleus) activities (Riechers etal. 2010, 2011; 
Danielson et al. 201 1). 

Owing to the negative K-correction, the number of submil- 
limeter galaxies (SMGs) does not need to peak at z~l-3, in- 
deed, if they exist, they should be as easy to detect at z=8 as 
they are at z=l-3. Their observed dearth is probably due to evo- 
lution of gas fraction and halo mass (Lacey etal. 2010; Hickox 
et al. 201 1). To be able to make a census of z~ 5 SMG and com- 
pare their number with lower-z SMGs, we use Herschel to iden- 
tify candidate high-z SMGs and strong lensing to provide the 
signal-to-noise and a unique spatial resolution. This technique 
has now proven quite successful in the study of high-z SMG 
(e.g. Smail et al 1997; Pello etal. 1999; Swinbank etal. 2010; 
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Lestrade etal. 2010; Negrello etal. 2010; Lupu etal. 2011; Cox 
etal. 2011). 

Strongly lensed sources are quite rare on the sky, and 
are serendipitously discovered in large surveys, such as the 
South Pole Telescope survey (Vieira etal. 2010), the Herschel 
HERMES (Oliver et al. 2010), or the huge area Herschel-ATLAS 
(H- ATLAS) survey (Eales etal. 2010). We have undertaken a 
more targeted survey with Herschel, focusing on galaxy clus- 
ters. Our key project, "Herschel Lensing Survey" (HLS-deep), 
has the goal to deeply image 44 massive clusters of galaxies with 
PACS and SPIRE (Egami et al. 2010; Rex et al. 2010), and con- 
tinues with a SPIRE snapshot survey (HLS-snapshot) of ~ 300 
clusters. The unique aim of the HLS survey is to find cluster- 
lensed sources, as opposed to galaxy-lensed objects found in 
field surveys such as H- ATLAS. These cluster-lensed sources 
suffer less from differential amplification than galaxy-lensed 
sources, while the dense environment of clusters increases the 
probability of galaxy lensing (Smail et al 2007). 

The HLS-deep survey is now fully executed and the source 
behind Abell 773 is the strongest source at 500/mi, among the 
44 clusters surveyed. Together with the z = 2.78 lensed galaxy 
in the Bullet cluster (Rex et al. 2010), it is one of only two HLS- 
deep SPIRE sources with a peak flux density at 500/im above 
100 mJy. It is located far from the cluster center (~4.5' = 1 Mpc), 
and likely lensed by a single galaxy. 

We present in this letter the continuum and molecular line 
observations of the Abell 773 SMG, deriving its gas excitation 
and physical properties. §2 reports the discovery of the SMG 
from our SPIRE and SMA data, IRAM spectroscopic observa- 
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tions to derive its redshift, with detections of CO, C i, H2O, [Nil] 
lines and CO(2-l) with EVLA. In §3 we model the data to infer 
the dust and gas content and its general properties. For distances 
we adopted the A-CDM concordance model: Ho = 71 km/s/Mpc, 
Q.m = 0.27 and D.\ — 0.73. 

2. Observations 

2.1. Herschel and SMA 

SPIRE observations were conducted using the "large-map 
mode", with orthogonal cross-linked scans producing 17'xl7' 
coverage in each of the three SPIRE bands. The maps include 20 
repetitions with a total observing time of 1 .7 hours, giving an on- 
source integration time of ~2000 s. The data were reduced using 
the Herschel interactive pipeline environment (HIPE v5.0). The 
maps were produced with the standard naive map maker, but in- 
clude additional turnaround data with scan speeds > 0.5"/s to 
increase the depth of the outer regions of the maps. For more 
details see Egami et al. (2010). 

HLSJ091828.6+514223 was discovered in the SPIRE map 
of the z=0.22 cluster Abell 773 (Barrena et al. 2007), shown 
in Figure Q] Its position is 0:2000 = 09 /! 18"'28.6 s and £2000 = 
51°42'23.3", too far away to be covered by PACS or Spitzer im- 
ages of the cluster. It lies close to an optical source (FiglTJ, which 
is likely the lensing galaxy. HLSJ09 1828.6+5 14223 has an SED 
(spectral energy distribution) peaking at 500/im, with F5oo=203 
mJy, suggesting a high redshift. The mm/submm contribution of 
the foreground galaxies is negligible. 

The Submillimeter Array (SMA) observed the source at 235 
GHz (1.3 mm) on 15 Jun 2011 with precipitable water vapor 
(pwv) between 1.5 and 4mm during 62min on-source, with 8 
GHz total bandwidth. The source was reobserved at 341 GHz 
(880/im) on 9 Dec 2011, during 330 minutes on-source and 
pwv of 1.5mm. Phase and amplitude gain variations with time 
were monitored using quasars 0927+390 and 0920+446, and the 
flux scale was determined using Titan (Table [TJ. The source is 
slightly resolved at 1.3mm by the 4.6x2.8" beam, and resolved 
at 880//mby the 2. 1x2.0" beam, with an N/S extension (FigQ}. 

2.2. IRAM and EVLA 

To measure the redshift of the source, we used the strategy de- 
veloped by Weiss et al. (2009) and Lestrade et al. (2010) with the 
multi-band heterodyne receiver EMIR at the IRAM 30-m tele- 
scope. At the time of observations (26-27 Sep and 6 Oct 2011), 
one setup with the 3mm receivers provided 7.43 GHz of instanta- 
neous, dual linear polarization bandwidth. We then covered the 
35 GHz, from 80 to 115 GHz, in five setups with some over- 
laps, integrating about one hour per setup, in wobbler-switching 
mode. The precipitable water vapor was between 2 and 5mm, 
and the system temperatures was 90-100K at 3mm and 200K at 
1.3mm. We used two backends, the 4MHz filter-banks, and the 
WILMA autocorrelator, which provides a spectral resolution of 
2 MHz. Pointing and focus offsets, determined once every two 
hours, were more accurate than 3". Continuum levels were ob- 
tained at 2mm and 1.3mm (and an upper limit at 3mm, cf Fig. 
|2]and TableQ]). At 1.3mm, we detected 40±5mJy, slightly lower 
than, but still compatible with the SMA value. 

The search revealed a line only at the fourth setup, at 92.3 
GHz. The fifth setup gave a second line at 1 10 GHz, suggesting 
a redshift, which was then confirmed by tuning the 2mm receiver 
at 129 GHz, the three detected lines thus being 12 CO J=5-4, J=6- 
5 and J=7-6 transitions, indicating a redshift of 5.2429. Other 
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Fig. 1. SPIRE 250/um (green, beam 18"), 500pm (red, 37") and 
SMA 880pm (black, 2") overplotted on the Subaru (suprime 
cam) V-image of HLSJ091828.6+514223. The latter shows the 
lensing galaxy at the center, slightly north of the SMA peak. The 
SMA contours start at 1 cr and increase in steps of 4<x. The EVLA 
CO(2-l) contours are peaked at the same position and have the 
same N-S elongation. 



lines were searched for, such as dense tracers HCN and HCO + 
at 3mm, water at 2mm, and the [NII]205pm line at 1.3mm. 

Around 1 hr of integration was also obtained in CO(2-l) us- 
ing the Expanded Very Large Array (Perley et al. 201 1), with ob- 
servations and data reduction following those reported by Ivison 
etal.(2011). 

Figure [3] displays all detected lines, and results of the 
Gaussian fits are reported in Table [2] 

3. Results 

3.1. Lens modeling 

The V and i images obtained with the Suprime Cam at Subaru 
shows the presence of a foreground galaxy at ~ 1 .2" north of 
the SMG source. A 3xl0min spectrum was taken on Keckll 
with DEIMOS 6001/mm grating at 7200A with the GG455 fil- 
ter, centered on the foreground galaxy using a 1" slit at PA=14° 
in ~1" seeing. The spectroscopic redshift of the lens was mea- 
sured as z=0.63±0.005, from Call-HK lines and the G band. 
Adopting an apparent size for the lensed source of ~ 2" from the 
SMA maps, we have used lenstool (Kneib etal. 1993) to pro- 
duce Einstein rings or crosses of this size, and estimate a mag- 
nification fi ~1 1 (sum of the images) for a lens at z=0.63. Since 
the images are not identified yet, we estimate that 1 1 is the max- 



Table 1. Flux S v (mJy) for HLSJ091828.6+514223, with IRAM- 
30m (I), SMA (S) and Herschel (H) 



3 mm 2mm 1.3mm 0.88mm 500/jm 350/mi 250//m 
I I S S H H H 

<2 15±7 55+7 125+8 203+9 168+8 85+8 



Not including absolute flux calibration uncertainties 
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Table 2. Observed line parameters toward HLSJ09 1828.6+5 14223 



Line 


^obs 


J- mb 


c 

J v 


AVfwhm 




V 


/'/in 10 








Lmts.J 


r™ T.rl 

LmJyj 


[kms 


- 1 1 

] 


[ Jy km s ^ ] 


L Kin s j 


[K km s ' pc^] 




CO(2-l) blue 


36.9280 




3.1 ± 0.7 


230 ± 


60 


0.8 ±0.1 


-530 ± 20 


17. ± 2. 


EVLA 


red 






4.0 ± 0.7 


510 ± 


30 


2.2 ± 0.1 


-20 ± 20 


55. ± 3. 




CO(5-4) blue 


92.3077 


1.7 ±0.3 


8.6 ± 1.5 


110 ± 


30 


1.0 ±0.2 


-540 ± 12 


4. ± 1. 


IRAM-30m 


red 




2.7 ± 0.3 


13.4 ± 1.5 


540 ± 


40 


7.7 ± 0.5 


-5+ 17 


30. ± 2. 




CO(6-5) blue 


110.7615 


1.7 ±0.4 


11. ±2 


160 ± 


60 


1.9 ±0.4 


-540 ± 25 


5. ± 1. 


IRAM-30m 


red 




3.1 ±0.4 


16 ± 2. 


510 ± 


50 


8.3 ± 0.6 


± 17 


23. ± 2. 




CO(7-6) blue 


129.2111 


2.0 ± 0.3 


9.7 ± 1.5 


150 ± 


20 


1.5 ± 0.2 


-510 ± 30 


3. ± .5 


IRAM-30m 


red 




2.9 ± 0.3 


14.3 ± 1.5 


560 ± 


20 


8.5 ± 0.3 


4 ± 10 


17. ± 1. 


and PdBI 


Ci( 3 P 2 -5Pi)blue 


129.6422 


1.6 ±0.3 


8.3 ± 1.5 


150 ± 


20 


1.4 ± 0.2 


-530 ± 9 


2.8 ± .5 


IRAM-30m 


red 




1.7 ± 0.3 


8.8 ± 1.5 


530 ± 


40 


4.9 ± 0.3 


20 ± 14 


10. ± 1. 


and PdBI 


H 2 O p (2,0,2-l,l,l) 


158.2481 


2.2 ± 0.5 


11 ± 2 


400 ± 


60 


4.8 ± 0.7 


-17 ±30 


6.5 ±0.8 


IRAM-30m 


Nil blue 


234.0469 


1.6 ± 0.2 


8 ± 1 


410 ± 


90 


3.5 ±0.6 


-620 ± 35 


2.2 ± 0.4 


IRAM-30m 


red 




1.7 ±0.2 


8.7 ± 1 


340 ± 


60 


3.2 ± 0.5 


-33 ± 30 


2.0 ± 0.4 





Data and errors are derived from Gaussian fits for the two V-components. 
* The velocity is given relative to z=5.2429 (centered on the red component of the CO(6-5) line) 
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Fig. 2. SPIRE (250, 350 and 500/rai), SMA (0.88 and 
1.3mm) and IRAM (1.3, 2 and 3mm) photometric data for 
HLSJ091828. 6+514223, superposed on our model (red line) of 
a single big grain component at T</=52K (J3 = 2). The M82 and 
Arp220 SED from NED are also plotted. Our source has an SED 
very similar to Arp 220. 



imum amplification, whose range is 5-11. The cluster does not 
boost the magnification by more than 10% (Richard et al. 2010). 

3.2. Dust emission 

Table Q] displays the available photometric data from SPIRE, 
SMA and the IRAM-30m. We fitted the data in Fig.f2]with a sin- 
gle optically thin big grain component at an average temperature 
7^=52K, an emissivity slope /3=2, and a mass absorption coeffi- 
cient k-50 cm 2 /g at 100/Kiri. We derived the total FIR(8-1000yum) 
luminosity Lfir = 1.1 10 l4 //j. L , and the star-formation rate SFR 
= 1.8 10 4 /^/ M Q /yr, applying the relation from Kennicutt (1998). 
Fitting with other dust models, such as Chary & Elbaz (2001), 
leads to the same luminosity, within the uncertainties. The mass 
of the dust is then Mj„ s/ = 5.0 10 9 /^ M , which for a gas-to-dust 
mass ratio of 150 leads to M gas = 7.6 10 n /j" M Q . Other fits are 
possible, although less good, with lower values of j3 (down to 
1 .5), higher dust temperatures and lower dust masses, down by 
a factor 4. 




10 I I I I I I i i — I I I I I I 10 I I I I I I I I I I I I I I I I 
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Fig. 3. Four CO lines, with the Ci( 3 P 2 - 3 Pi), H 2 O p (2,0,2- 1,1,1) 
and [Nil] lines. All spectra are from the IRAM-30m, except for 
the CI/CO(7-6) from the PdBI, and CO(2-l) from EVLA. Upper 
vertical lines indicate the red (R) or blue (B) components for 
CI, and lower vertical lines indicate the same for all others. All 
panels have the same velocity scale (translated for CO(5-4) and 
CO/CI). The blue component of the [Nil] line is much stronger 
and broader than for the CO lines, suggesting a strong ionized 
gas flow. Exact velocities are detailed in Table [2] 
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Fig. 4. Peak CO line fluxes for the red (dots) and blue compo- 
nent (triangles), with the best LVG model (bold and light dashed 
lines, respectively) computed with a density nn 2 =3.5 10 3 cm -3 , 
T* = 45K, and a column density N(CO) = 10 18 cm^/kms" 1 . 
The light dashed line is related to the bold line in a ratio of 0.8 
to 1. The color lines represent CO data from other galaxies (e.g. 
Weiss etal.2007; Danielson et al. 2011). 



3.3. CO lines 

Four CO lines are clearly detected in two components, as shown 
in Figure[3] The peak flux of the stronger red component is mod- 
eled in Figure [4] with an LVG (large velocity gradient) code, 
varying the H2 volume density, the CO column density and the 
kinetic temperature (Combes etal. 1999). The blue component 
has a peak intensity lower by 25%, whatever the J level, within 
the error bars. The three high-J CO lines are well fitted by a sin- 
gle excitation component, peaking at the 7 = 6 level. We found 
the best fit with a density nn 2 = 3.5 10 3 cirT 3 , T* = 45K, and 
a column density N(CO) = 10 18 crrr 2 /kms _1 . The solution is 
somewhat degenerate, between low H2 density and higher CO 
column density, and we favor optically thick CO emission. PDR 
(photodissociation region) models fit as well, the CO being emit- 
ted by small and dense optically thick clouds. Observations of 
higher J CO lines are required to distinguish between PDR and 
XDR (Xray-dissociated region), where CO is emitted by a more 
extended and diffuse hot medium (Meijerink etal.2007). The 
observed CO(2-l) flux is 50% higher than predicted by the high- 
excitation component, suggesting the presence of a cooler com- 
ponent, containing one third of the mass. Adopting the conver- 
sion factor M(H 2 )/L^ O(1 _ 0) =0.8 M /(K kms^pc 2 ) of ULIRGs 
(ultraluminous galaxies, e.g. Solomon etal. 1997), we derive 
from the CO(2-l) integrated flux an H2 mass of 5.8 I0 n /fi 
M Q . This mass is 0.8 times that derived from the dust emission, 
but there are still uncertainties in the conversion ratio to apply. 
The derived ratio Lfir/L^ 0(3 _ 2) ~ 260, compatible with other 
(U)LIRGs (Iono et al. 2009). 

3.4. Atomic carbon 

Since the C i( 3 />i - 3 Pq) line is not available in the receiver range, 
we have observed the Ci^i^-^i) with the 30m, and confirmed 
it later with the Plateau de Bure (PdBI), with a better signal- 
to-noise ratio (Boone etal. in prep). The latter is shown in Fig 
[3j with the two velocity components detailed in Table [2] From 



an adopted excitation temperature of T ex — 40K, and a total lu- 
minosity L' a(3p _ 3p j= 12.8 10 10 K kms~'pc 2 , we derive a car- 
bon mass of M a =2 10 8 //i M (Weiss etal. 2005). This yields 
a C1/H2 number abundance of 6+2 10~ 5 , comparable to other 
high-z SMGs (e.g. Danielson etal. 201 1). 



3.5. Water 

ISO and Herschel observations have revealed that far-infrared 
H2O lines are surprisingly strong with respect to CO in ULIRGs, 
in particular those containing an AGN (Arp 220, Cloverleaf, Mrk 
231, Gonzalez-Alfonso etal. 2004, 2008; Bradford etal. 2009; 
van der Werf etal. 2010). Water is thought to evaporate from 
grains in shocked regions, in dense hot cores due to cosmic rays 
or possibly X-rays near AGN (Gonzalez-Alfonso et al. 2010). In 
Mrk231, van der Werf etal. (2010) have shown that starbursts 
and PDR alone cannot explain the high-J excitation of the CO 
lines, nor the H2O/CO ratios, and a very strong FIR radiation 
field or XDR conditions related to AGN are required. The asso- 
ciations of starbursts with AGN might be even more frequent at 
high redshift, and a strong H2O line has been detected at z=2.3 
(Omont etal. 2011) and z=3.91 (e.g. van der Werf et al. 201 1). 

The water line detected in HLSJ09 1828.6+5 14223 is the 
highest flux line expected from the mixed model of Gonzalez- 
Alfonso & Cernicharo (1999), including collisional excitation 
and far-infrared pumping. We did not detect the H-20 ( 1,1,0- 
1,0,1), with an upper limit 3 times lower than the H2O p (2,0,2- 
1,1,1) line, which is compatible with the model (expected ra- 
tio of 4). The water lines are very optically thick, but infor- 
mation comes from their ratio to the CO lines, which here is 
I(H 2 0)/I(CO) ~ 0.5, equal to the ratio observed in Mrk231 or 
the Cloverleaf quasar at z=2.6, SDP 17b at z=2.3, and a little 
lower than in APM08279 at z=3.91 (van der Werf etal. 2011, 
and references therein). In all these galaxies, this ratio is 1 - 
2 orders of magnitude higher than in Orion bar. This sug- 
gests the existence of a hot dense component, or an XDR in 
HLSJ091828.6+514223. 



3.6. The [Nll]205fim line 

Figure[3]displays the spectrum of the 3 Pi-> 3 Po ground state fine 
structure line of ionized nitrogen, smoothed to a resolution of 
100 km s -1 . The two velocity components are clearly detected, 
and the blue component is much broader than for the CO and 
CI lines. This fine structure line, with a low excitation poten- 
tial, mainly traces HII regions photoionized by OB stars. In the 
Milky Way, the line appears moderate, with an intensity ratio 
[CII]/[NII] ~ 10. The latter ratio is 20 in M82 (Petuchowski 
etal. 1994), and 15 in Mrk231 (Fischer etal. 2010). Only up- 
per limits have been reported in high-z galaxies, the most con- 
straining being in the Jl 148+5251 quasar at z=6.4 (Walter 
etal. 2009), with L [NH] /L FIR < 2 10~ 5 , and L [Nn] fL CO (6-5) < 
0.9. Here we detect the line with L^/j/L/t/r = 4 10~ 5 and 
L[tv//]/Lco(6-5) = 0.15. These values are intermediate between 
the AGN-dominated Mrk231, where L [NU] fL F!R = 0.8 10~ 5 , and 
the starburst M82, where L^jj/Lfm = 7.5 10~ 5 . The broad 
line seen in the [Nil] blue velocity component could be re- 
lated to galactic winds, often seen in AGN/starburst compos- 
ite ULIRGs, and with higher velocities than in neutral flows 
(Veilleux et al. 2005; Rupke et al. 2005; Fischer et al. 2010). 
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Table 3. Line upper limits at 3<x in channels of 140 km s 



Line 


v obs (GHz) 


S v (mJy)* 


HCN&HCO+(6-5) 


-86. 


< 4. 


13 CO(5^1) 


88.2484 


< 10. 


H 2 (1, L0- 1,0,1) 


89.2111 


< 4. 


CS(ll-lO) 


86.2882 


< 10. 



3.7. Other lines 

We have searched for the high-density tracers HCN and HCO + , 
in their lowest available level, i.e. J=6-5 (Table [3). The upper 
limits found HCN(6-5)/CO(6-5) < 0.25 are not yet constraining, 
since ULIRGs have on average HCN/CO ~ 0.1, except for 1 or 
2 luminous quasars (Gao & Solomon 2004). 

4. Discussion and conclusion 

We have estimated an amplification //ill from lenstool (likely 
range 5-11). The amplification factor can also be estimated 
from the trend of the CO(1-0) FWHM line width with lu- 
minosity L (Harris etal. 2012; Bothwell etal. 2012 in prep): 
(j, = 3.5*(L'/10 11 ^kms" 1 pc 2 )*(^W//M/400kms" 1 r L7 . This 
yields ju=16±5, which is consistent with our earlier estimate, and 
suggests that the amplification is likely ~ 1 1. Taking into account 
this magnification, the intrinsic FIR luminosity of the source is 
still in the ~ 10 13 L G range, corresponding to a hyper-luminous 
source. The detected lines (except water) are divided into two ve- 
locity components, which suggests a galaxy merger, and would 
explain the high derived SFR (Engel et al. 2010). 

At these high luminosities, an AGN component is probable, 
which is supported by the high water emission. However, the 
flux in the CO ladder peaks at J - 6 as in starbursts like NGC 
253 and M82, and we have no evidence yet that it is extended to 
higher J as in powerful AGNs like Mrk231 or APM08279. The 
[NII]205/im line is detected for the first time at high-z, and its 
relative luminosity with respect to Lf /fl is intermediate between 
a starburst (M82) and an AGN (Mrk231). Its blue velocity com- 
ponent is much broader than for the CO lines, which could come 
from violent flows of ionized gas. The two velocity components 
have quite different properties (in H2O, [Nil]), favoring an inter- 
pretation in terms of two galaxies, with a high velocity difference 
of 540 km s . Observations of the higher-J CO lines are needed 
to complete the gas excitation, and interferometer maps of key 
lines, such as the [CII] line, will bring more information about 
the nature of this interesting system. 
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